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We compute and analyze the Flavor-Changing Neutral Current (FCNC) interactions of Minimal Supersym- 
metric Standard Model Higgs bosons (h = h°, H°, A ) with heavy quarks (top and bottom), focusing on 
the strongly-interacting sector. We correlate the Higgs bosons production cross-section at the LHC with the 
FCNC decay branching ratios and find the maximum allowed values of the production rates, cr(pp — > h — > 
qq') = o{pp — ► h) x B(h — > qq') (qq' = tc or 6s) after taking into account limits from low energy data 
on flavor-changing interactions. We single out the top quark channel, with a maximum production rate of 
<r max (pp — ► h — ► tc) ~ 1CT 3 — 10" 2 pb, as the most promising FCNC channel to be detected at the LHC. 



1. Introduction 

The possibility that Nature hides the last pos- 
sible symmetry of the S'-matrix, Supersymmetry 
(SUSY), in the form of a broken symmetry of the 
fundamental interactions is one of the most in- 
triguing puzzles over the last 30 years of high 
energy physics. SUSY offers a most tantalizing 
paradigm to unify Particle Physics interactions 
with Gravity, and it might even provide a clue 
to a more fundamental superstring unification of 
all interactions. As such SUSY will be scruti- 
nized in great detail at the LHC. It also offers 
a possible solution to the longstanding natural- 
ness problem in the Higgs sector of the SM [1]. 
If SUSY is realized around the TeV scale, the 
LHC experiments shall be able to directly pro- 
duce the SUSY particles for masses smaller than 
a few TeV [2]. Among the many strategies de- 
vised to hunt for SUSY, Flavor- Changing Neutral 
Currents (FCNC) offer perhaps a unique labora- 
tory to seek for the new signs of physics beyond 
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the SM. In the SM the FCNC effects are com- 
pletely absent at the tree-level. At one loop, how- 
ever, they are possible, but then the contributions 
from the new particles enter on equal footing with 
those from SM particles. In certain regions of the 
parameter space the new physical effects may well 
dominate the SM contributions. This is particu- 
larly so when the SM one-loop effects turn out 
to be highly suppressed. In such situations the 
sole observation of these FCNC processes would 
be instant evidence of new physics. An example 
of this kind of appealing scenarios occurs within 
the FCNC physics of the SM Higgs boson (Hsm) 
interactions with quarks. For example, the FCNC 
vertex Ho,yitc may lead (at one loop) to such rare 
decays as t — ► Hsm c or Hsm — > t c, depending on 
the mass of Hsm- Both of these modes are ex- 
tremely suppressed at one loop, with branching 
ratios of order 10~ 14 or less [3,4], hence 10 or- 
ders of magnitude below other more conventional 
(and relatively well measured) FCNC processes 
like b — > sj [5] . 

The Minimal Supersymmetric Standard Model 
(MSSM) introduces new sources of FCNC inter- 
actions mediated by the strongly-interacting sec- 
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tor 1 . They are produced by the misalignment 
of the quark mass matrix with the squark mass 
matrix, and the main parameter characterizing 
these interactions is the non-flavor-diagonal term 
in the squark-mass-matrix, which we parame- 
terize in the standard fashion [8] as (M )y = 
5ijfnifhj (i 7^ j), rhi being the flavor-diagonal 
mass-term of the z-flavor squark. Since there are 
squarks of different chiralities, there are different 
Sij parameters for the different chirality mixings. 
In this work we will assume flavor-mixing only 
among the Left-Chiral squarks, since these mix- 
ing terms are expected to be the largest ones by 
Rcnormalization Group analysis [9]. 

Some work in relation with the MSSM Higgs 
bosons FCNCs has already been performed [4, 7, 
10,11]. In this work we compute and analyze 
the production of any MSSM Higgs boson (h = 
h°,H°, A ) at the LHC, followed by the one-loop 
FCNC decay h — > b s or h — > tc, and wc find 
the maximum production rates of the combined 
cross-section: 



cr(pp -> h -> qq') 



B{h 



qq) = 



a(pp -> hX)B(h -^qq' 
qq' + qq') 



(1) 



bs 



qq' being a pair of heavy quarks (qq' 
tc), taking into account the restrictions from the 
experimental determination of B (b — ► 57) [5]. For 
other signals of SUSY FCNC at the LHC, without 
Higgs bosons couplings see Ref. [6] . 

2. Computation setup 

In this section we give a summarized explana- 
tion of the computation. For further details see 
Refs. [4,11]. 

We include the full one-loop SUSY-QCD con- 
tributions to the FCNC partial decay widths 
T(h — ► qq') in the observable 0). 

The Higgs sector parameters (masses and CP- 
even mixing angle a) have been treated using the 
leading m t and rrib tan approximation to the 
one-loop result [12]. 

The Higgs bosons total decay widths T(h — + X) 
are computed at leading order, including all the 

1 For a brief description of these interactions see e.g. 
Ref. [6], for details see [7]. 



relevant channels: T(h -> ff,ZZ,W + W~,gg). 
The off-shell decays T(h -> ZZ* ,W ± W T *) have 
also been included. This is necessary to consis- 
tently compute the total decay width of T(h° — > 
X) in regions of the parameter space where the 
maximization of the cross-section Q is obtained 
at the expense of greatly diminishing the partial 
decay widths of the two-body process h° — > bb 
(due to dramatic quantum effects that may re- 
duce the CP-even mixing angle a to small val- 
ues [13]). The one-loop decay rate T(h — ► gg) 
has been taken from [14] and the off-shell decay 
partial widths have been recomputed explicitly. 

The MSSM Higgs boson production cross- 
sections have been computed using the programs 
HIGLU 2.101 and PPHTT 1.1 [14,15]. These 
programs include the following channels: gluon- 
gluon fusion, and associated production with top- 
and bottom-quarks. We have used the leading 
order approximation for all channels. The QCD 
renormalization scale is set to the default value for 
each program. We have used the set of CTEQ4L 
PDF [16]. 

For the constraints on the FCNC parameters, 
we use B(b -> 37) = (2.1 - 4.5) x 10~ 4 as the 
experimentally allowed range within three stan- 
dard deviations [5] . We also require that the sign 
of the b — > sj amplitude is the same as in the 
SM [17] 2 . 

Running quark masses (m q (Q)) and strong 
coupling constants (a s (Q)) are used throughout, 
with the renormalization scale set to the decaying 
Higgs boson mass in the decay processes. 

Given this setup, we have performed a Monte- 
Carlo maximization [18] of the cross-section JJJ 
over the MSSM parameter space, keeping the pa- 
rameter tan P fixed, and under the simplification 
that the squark and gluino soft-SUSY-breaking 
parameter masses are at the same scale: 



m~, = raj = m,~. 



u a — m g = M3USY 



(2) 



Although we have performed a full one-loop 
computation, involving the exact diagonalization 
of the 6x6 squark mass matrix, it is enlightening 
to look at the approximate leading expressions to 



2 This constraint automatically excludes the fine-tuned re- 
gions of Ref. [4] . 
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Table 1 

Maximum values of B{h — > bs) and corresponding 
SUSY parameters for m^o = 200 GeV, tan j3 = 50. 



Figure 1. The maximum value of B{h — > bs) as a 
function of m^o, for tan (3 = 50. 



understand the qualitative trend of the results. 
The SUSY-QCD contribution to the b -> S7 am- 
plitude can be approximated to: 

A s ^(& 1 S7 ) ~ S 23 m ^- A ^P) f (3) 

-^SUSY 

whereas the MSSM Higgs bosons FCNC effective 
couplings behave as: 

f sin(/3-a cff ) (if ) 

9h q ?~S 23 -JP^{ cos(/3-a eff ) (ft ) .(4) 
m susy I 1 (^0-) 

The different structure of the amplitudes in 
eqs. (|3J| and will permit to obtain and ap- 
preciable FCNC Higgs boson decay rate, while 
the prediction for B(b — > sj) stays inside the ex- 
perimentally allowed range. 

3. Bottom-strange channel 

For the analysis of the bottom-strange pro- 
duction channel, we study first the Higgs boson 
branching ratio . Figure 2] shows the maxi- 
mum value of B{h — * bs) as a function of the 
pseudoscalar Higgs boson mass m^o . We observe 
that fairly large values of B(h° — > bs) ~ 0.3% are 
obtained. Table 2] shows the actual values of the 
maximum branching ratios, and the parameters 
that provide them for each Higgs boson. Let us 



discuss first the general trend, which is valid for 
all studied processes: the maximum is attained at 
large M SUSY and moderate <5 23 . The SUSY-QCD 
contribution to b — ► 57 (|3J) decreases with -MguSYj 
therefore to keep B(b — > sj) in the allowed range 
when Msusy is small, it has to be compensated 
with a low value of 623, providing a small FCNC 
effective coupling Q. On the other hand, at large 
Msusy the second factor in eq. J3J decreases, al- 
lowing a larger value of £23- Thus, the first factor 
in eq. (@J grows, but the second factor in eq. J3J 
stays fixed (provided that ~ A^susy), overall 
providing a larger value of the effective coupling. 
On the other hand, a too large value of #23 has 
to be compensated by a small value of ^l/^SUSY 
in provoking a reduction in ijlf. In the end, 
the balance of the various interactions involved 
produces the results of Table ^ 

The maximum value of the branching ratio for 
the lightest Higgs boson channel is obtained in 
the small a e s scenario [13]. In this scenario the 
coupling of bottom quarks to h° is extremely sup- 
pressed. The large value of B(h° — > bs) is ob- 
tained because the total decay width T(h° — » V) 
in the denominator of 0) tends to zero (Ta- 
ble Q), and not because a large FCNC partial de- 
cay width in its numerator [4]. 

The leading production channel of h° at the 
LHC at high tan/3 is the associated production 
with bottom quarks, and therefore the h pro- 
duction will be suppressed when B(h° — > bs) is 
enhanced. We have to perform a combined analy- 
sis of the full process to obtain the maximum 
production rate of FCNC Higgs bosons meditated 
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Figure 2. Maximum SUSY-QCD contributions to 
cr(pp — > h — > b s) as a function of m^o , tan (3 = 50. 



Table 2 

Maximum value of <r{jpp — > /i — > 6 s) at the LHC, 
for m^o = 200 GeV and tan/3 = 50. Shown are 
also the corresponding values of B(h — > 6s) and 
of the total width of the Higgs bosons, together 
with the values of the SUSY parameters. 



events at the LHC. Figure [5] and Table |2] show 
the result of the maximization of the production 
cross-section (QJ. The central column of Table 
shows that when performing the combined maxi- 
mization T(h — » X) has a much larger value, and 
therefore the maximum of the combined cross- 
section is not obtained in the small a e ff scenario. 
The number of expected events at the LHC is 
around 50,000 events/ 100 fb" 1 . While it is a large 
number, the huge 6-quark background at the LHC 
will most likely prevent its detection. Note, how- 
ever, that the maximum FCNC branching ratios 
are around 10 _4 -10 -3 , which is at the same level 
that the already measured B(b — > sj). 

4. Top-charm channel 

The results of the numerical scan for this chan- 
nel are similar to the 6s channel, so we will focus 
mainly on the differences. Figure |31 shows the 
maximum value of the production cross-section 
o~(pp — > h — > tc) as a function of m^o, while 
Table [3] shows the actual value of the maximum, 
together with the SUSY parameters that provide 
them, for m^o = 300 GeV. Only the heavy neu- 
tral Higgs bosons contribute to this channel. The 
general trend explained in section [3] is also valid 



here 3 . From Table 01 we would expect some 300 
events/100fb _1 at the LHC. However, the maxi- 
mum is attained at low tan/3, and we have fixed 
a moderate value of tan/3 = 5. By taking lower 
values of tan /3 we find that the number of events 
grows fast [11], e.g. up to - (500,900,2000) if we 
would have chosen tan j3 = (4, 3, 2) respectively. 
Due to the single top quark signature the detec- 
tion of this channel should be feasible at the LHC. 

5. Conclusions 

We have performed a computation of Higgs bo- 
son mediated FCNC events at the LHC. Alterna- 
tive effects of SUSY FCNC have been considered 
in Ref. [6]. Note from eq. (0} that the present 
effects, though smaller than the latter, do not 
decouple with Msusy = frig = |/i|. The maxi- 
mum expected branching ratio for the 6s channel 
is large, and we expect a maximum of ~ 50, 000 
events/100 fb -1 in this channel at the LHC, which 
might be difficult to detect due to the huge back- 
ground. For the tc channel we expect a maximum 
of cr max (pp -> h -> tc) ~ 10~ 3 - 10~ 2 pb, which 
means several thousand events per 100 ftT 1 at the 

3 Recall that the Sij parameters in the up-sector are related 
to the corresponding parameters in the down-sector by the 
Cabibbo-Kobayashi-Maskawa matrix, see e.g. [8], and are 
therefore constrained by B(b — » 57). 
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Figure 3. Maximum SUSY-QCD contributions to 
o~(pp —> h — > tc) as a function of m^o , tan (3=5. 



LHC. Due to the single top quark signature they 
should be easier to detect than the bs channel, 
providing the key to a new door to study physics 
beyond the Standard Model. It is now an experi- 
mental challenge to use this key to open the door, 
and prove that these events can be effectively be 
separated from the background. 
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